Asymmetric coherent transmission for single particle diode and gyroscope 
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We study the single particle scattering process in a coherent multi-site system consisting of a tight-binding 
ring threaded by an Aharonov-Bohm flux and several attaching leads. The asymmetric behavior of scattering 
matrix is discovered analytically in the framework of both Bethe Ansatz and Green's function formalism. It is 
found that, under certain conditions, a three-site electronic system can behave analogous to a perfect semicon- 
ductor diode where current flows only in one direction. The general result is also valid for a neutral particle 
system since the effective magnetic flux may be implemented by a globe rotation. This observation means that 
the three-site system can serve as an orientation measuring gyroscope due to the approximate linear dependence 
of the current difference of two output leads on the rotational angular velocity. 
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I. INTRODUCTION 

Single particle quantum devices work essentially in a quan- 
tum mechanical way by using the whole features of quantum 
states, especially the phases of quantum states. A typical ex- 
ample is the single electron transistor (SET) |H 0> HI H Htl 
and its photon analogue which plays a central role in 

quantum manipulation and quantum measurement. A SET is 
a mesoscopic system that allows confined electrons to tunnel 
to the metallic leads. It turns on and off again every time one 
electron is added to the isolated region ]2[]. Unlike the conven- 
tional transistors, which can be understood by using classical 
concepts, the SET is a quantum mechanical one in substance. 
It can be utilized for measuring the quantum effects in Joseph- 
son junction superconductor circuit and nano-mechanics res- 
onators 

In this article, we will pay attention to another single parti- 
cle quantum device, which can be understood as the quantum 
analogue of the conventional diode device. In practice, diode 
devices are indispensable for building various electronic cir- 
cuits in the traditional electronics. Accordingly, once a co- 
herent diode device has been implemented, a single-particle 
quantum circuit may be further realized in the level of single 
quantum state by using of the two basic elements, single par- 
ticle transistor and diode. Compared with the traditional elec- 
trocircuits based on the density distribution of electrons, the 
new quantum circuit makes full use of the quantum properties 
including quantum phases in information processing. In this 
sense, we can also regard it as a necessary element in quantum 
information science and technology. 

A major research effort in recent years is to seek a sim ples t 
single particle system with diode features 11 ll Il2l Il3l Il4ll . 
which is characterized by the asymmetric performance of 
transmission coefficients along the opposite directions. In the 
present investigation, we find that a three-site tight-binding 
ring system threaded by an Aharonov-Bohm (AB) magnetic 
flux has analogy to a perfect diode in a wide range of param- 
eters. 

In addition, our studies are not restricted to the charged par- 
ticle case. For a neutral particle system, an effective magnetic 



flux is induced by applying a globe rotation [15], which can 
also break the time reversal symmetry and lead to the asym- 
metric transmissions. It is shown that the difference of the two 
output currents is linearly proportional to the rotational angu- 
lar velocity approximately, which behaves as a new kind of 
extended Sagnac gyroscope lfl6l[T7ll . 

Decades ago, Biittiker et al. have found that the conduc- 
tance of a multi-terminal sample is asymmetric in the pres- 
ence of an AB flux 118l Il9l 12011 . Realization of such mech- 
anism in quantum device requires its size being smaller than 
the phase-breaking length. Thus seeking a minimized system 
with asymmetric feature is attractive in practice. A discrete 
system may be a good candidate to accomplish such a task. 
It is of both theoretical and practical importance to study this 
theory in a discrete system. Moreover, we show practical ap- 
plications of the quantum interference effect as single electron 
diode and single particle gyroscope. In the general discus- 
sion part, it is interesting to find that some special configura- 
tions may protect the transmission coefficients from symmetry 
breaking even if the AB magnetic flux is imposed. 

This paper is organized as follows. In Sec. II, the general 
Af-site ring-shaped model with attaching leads is presented, 
and an analytical method for calculating the transmission co- 
efficients is introduced. In Sec. Ill, we study the asymmetric 
transmission behavior of the 3-site system by using the ex- 
act solution. In Sec. IV, a single particle gyroscope and its 
physical realization in optical lattice is proposed. In Sec. V, 
we give a general discussion about the condition of asymmet- 
ric transmission by making use of Green's function approach, 
which also confirms the results obtained above. Conclusions 
are summarized at the end of the paper. In the appendix, 
a proof is given to show the equivalence between the Bethe 
Ansatz and the Green's function method. 



II. MODEL AND ITS EXACT SOLUTION 

The central system we concern is described by an A^-site 
tight-binding ring threaded by a magnetic flux shown in Fig. 
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FIG. 1: (color online) (a) Configuration of the ring shaped scattering 
system including M arc chains {a [1 ',a , . . . , a [M1 ) and M attaching 
leads {b ll \b l2 \ . . . , fo [M1 j threaded by a magnetic flux. The input and 
possible output currents are marked by wavy arrows, (b) Schematic 
illustration of the three-site chiral coherent scattering system with 
diode features. 



[TJa). The Hamiltonian reads 



He = -J / t {<■"'' <i ,<i i- : + H.c.) + 0) V fltay. 



(1) 



7=1 



7=1 



Here, aj (ay) is the fermion creation (annihilation) operator at 
the j'th site. The external magnetic field does not exert force 
on the Bloch electrons, but makes hopping integral J between 
sites j and (J + 1) pick up an AB phase factor exp (i(/>jj in the 
Peierls approximation iTJlll . where 



<Pj 



2n r i+l 

00 J j 



A-dl, 



(2) 



<po = hc/e is the flux quanta and <t> = Y/j =l 4>j is me t° ta l 
magnetic flux; a> is the chemical potential of the central sys- 
tem. In the same figure, M half-infinite tight-binding leads are 
attached to M sites on the ring, and the corresponding Hamil- 
tonian are 



Hr 



z*~z 

l=\ l=\ 



H.c. 



(3) 



Without loss of generality, we let I = l„ lead as the input lead, 
while all the other indexed / + l m leads are output ones. The 
spin degree of freedom is omitted for notational brevity here 
since the model does not contain interaction involving spin. 

In order to investigate the scattering problem through such 
a system in the framework of Bethe Ansatz |0, HH S] more 
efficiently, we regard the /V-site ring as M arc tight-binding 
chains connected by the head and tail. As illustrated in Fig. 
Q2a), the length of the /th arc chain a® is Ni, and the Zth lead 
is joined with the 1th site of the /th arc chain. The total 
number of sites on the ring satisfies N i = N - Then the 
Hamiltonian H = He + Hi is rearranged with 



M 



Hr = 



Hr = 



Ni-l 



1=1 [j=\ 

Ml \ 

-£ g a™bf + r£ b ™ b n + H.c. . 



(4) 



The scattering wave function is supposed to be 

tf? (/) = A, (0 W + A 2 (/) e <**Z& ^ ), 

ififU) = B(l)e ikJ + 6 kki e- ikj , (5) 
where B(lj„) actually defines the reflection coefficient 

RL. = mk,t, (6) 



while B (I) for I + k n give the transmission coefficient from 
the /„,th lead to the /th lead, 



ib cor 



(7) 



The coefficients A\ (I), A 2 (I) and B (I) should satisfy the con- 
necting conditions and the Schrodinger equation 

^- 1] (/V M +1)=^(1), 

e '<vw(2 )+e -<",vr i] (^-i)" 



= (E-w)^(l), 

-t>$(2)-gtf(l) = E^(l). (8) 

The relation between k and q is also given by the Schrodinger 
equation, 



-/^>W(; + i) + e -^,^ ] (;-i) 

-t[^(j+l) + ^(j-l)] = E^(j), 
for y # 1, i.e., 

E = -2t cos k = -27 cos q + a> 



= (£-«) ^(j), 



(9) 



(10) 



Then the 3M coefficients {Ai (0 , A 2 (0 , B (/)} are fully deter- 
mined by solving the above 3M independent equations in Eq. 
©. It shows that the Bethe Ansatz method can provide the 
exact wave function of the scattering state. For small size sys- 
tem, the explicit wave function allows us to clarify the mech- 
anism of the asymmetric transmission from the viewpoint of 
interference. 



IH. ASYMMETRIC COHERENT TRANSMISSION 

Hereafter, we mainly consider the simplest case of /V = 3 
and M = 3 shown in Fig. \Hb). The Bethe Ansatz approach 
gives the exact solution of transmission and reflection coeffi- 
cients as 



7*13 = 732 = T%\ = T R , 

T\2 - 723 = ?31 = T L , 

R\ — Ri — R3 = 1 - 7*12 _ 7^i3 = R, 



(11) 



where 



T R = 



^psin 2 fc|/r e - /o + 0|\ 



= 2t z cos k - g z e' k + tco, 

E = I2./Y cos <D + 37 2 f 2 © - © 3 



(12) 
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FIG. 2: (color online) Typical cases of transmission coefficients in 
the clockwise direction T R (black solid lines) and anti-clockwise di- 
rection T L (red dashed lines), (a) Tr and T L with respect to E/t in 
the case of t = J (upper panel) and t = J/2 (lower panel), where 
g 2 /t = J, O = 7r/2, and a) = 0. It shows that the system behaves 
like a perfect diode at E = (upper panel) and within a wide range 
E e [-/, t] (lower panel), (b) Transmission spectrum in the case of 
g = t = J and a> = with O = n/3 (upper panel) and (t> = 2n/3 
(lower panel). The perfect diode energy can be adjusted by the mag- 
netic flux O. (c) Plots of T R and T L at E = as a function of t/J with 
parameters g = t, O = n/2 and a> = 0. The transmission coefficients 
of two directions coincide in the limit of both t — > and t — > oo, but 
diverge dramatically at the critical point t = J. (d) When g 2 /t = J, 
i> = n/2, and E = 0, the diode feature is optimal at the resonant 
point o)j J = 0. 



We notice that in general, the transmission coefficients Tr 
and Ti from one lead to the other two are not identical, i.e., 
the current flow across two arbitrary leads is unidirectional. 
On the other hand, the clockwise coefficients T\j, 732, and T%i, 
defined as Tr, have the same value, so do the anti-clockwise 
coefficients T\2, T23, T31 defined as Tl- Namely, the transmis- 
sion coefficients have dextrorotary and levorotatory charac- 
teristics. We will discuss the general asymmetric scattering 
problem in detail in Sec. V. We will see that this kind of fea- 
ture is reasonable for a non-bipartite system with broken time- 
reversal symmetry induced by the external magnetic flux. 

Now we concentrate on some special cases of the 3-site ring 
to exemplify this feature. 

(1) g 2 It — J, O = n/2, and u> — 0, the Bloch electrons are 
injected at E — with momentum k = n/2. Eq. ( fTTT i and (fT2] i 
becomes 



IT, = 0,R = 0. 



(13) 



The profiles of the corresponding transmission spectra, Tr and 
Tl as functions of E for t = J and t = J/2, are plotted in Fig. 
|2]a). One can see that, when t = J/2 and g 2 /t = J, the diode 
is perfect even if the input energy E is shifted within the re- 
gion [—t, t]. This shows the great tolerant as a perfect quantum 
device. 

(2) g — t = J and 10 — 0, the input energy is determined by 
k as E — -2f cosfc, 



T R = l,T L = 0,R = 0, fork=7r-® + 2nn, 
T R = Q,T L = 1,R = 0, for it = <D -n + 2nn, 
T R = 0, T L = Q,R = 1, for £ = nn. 



(14) 



It means that the current can only flow in the clockwise direc- 
tion with no reflection. Thus the system is a perfect diode if 
we regard the 1th lead as the source and 2th lead as the drain. 



Therefore, for a given E, we can get a perfect diode device 
by tuning <£. The transmission spectrum of the <I> = n/3 and 
27T/3 cases are shown in|2jb) to give a visual impression. 

(3) g = t,<b = n/2 and u = 0, we plot T R and T L at E = 
with respect to t in Fig. [2jc). The point of t — J denotes the 
diode case. Beyond that critical point, the difference between 
Tr and Tl gets smaller. In the limit of t — > 0, Tr and Tl 
converged to 4/9. 

(4) g 2 /t = J, O = n/2, and E = 0. The dependence of T R 
and Tl on the frequency detuning a> is shown in Fig. |2jd). It 
is apparent that the resonant case is advantageous to forming 
a perfect diode. 

A diode system based on the above mechanism is differ- 
ent from the classical ones in semiconductor electronics. The 
conventional diode such as P-N junction is made up based on 
the different density distribution of electrons in the p-type and 
n-type materials. When the external voltage is absent and the 
equilibrium is reached, there is a difference of chemical po- 
tentials between the two materials. Consequently, current will 
flow readily in the forward biased direction since the applied 
voltage decreases the barrier; but not in the reverse biased di- 
rection because the barrier is raised. By contrast, the single 
particle diode device presented here works in the region with 
no chemical potential difference. It makes full use of the pure 
quantum interference phenomenon induced by the Aharonov- 
Bohm effect. 



IV. SINGLE PARTICLE GYROSCOPE 

So far we have seen that the threaded magnetic flux plays 
an important role in controlling the amount and direction of 
currents. It can be applicable to the more extended system. 
Actually, if the system is rotated, an effective magnetic field 
will be induced in the rotating frame of references. Therefore, 
the difference between the currents of the two output leads 
can reflect the rotational angular velocity of the system. That 
is the basic idea of making up a gyroscope by using the above 
mentioned asymmetric scattering system. The analysis below 
can be applied to a neutral-particle system by simply choosing 
= 0. 

For a rotating system with angular frequency O, an addi- 
tional term 



4 



N 



3.7+1 



+ H.c), (15) 



should be added on the Hamiltonian lfl5tl in the non-inertial 
frame, where K is a constant depends on the geometry of the 
central system. In the meantime, the Hamiltonian of leads Hi 
remains the same since the leads are along the radial direc- 
tions. Then the Hamiltonian of the central system becomes 
H c + H R or 



H' 



a.cij + \ 



+ H.c. 



(16) 



./=! 



where J n = V/ 2 + Q. 2 K 2 , tan^n = SIK/J. The total effective 
magnetic flux 



On = + 30n = # + 30n 



(17) 



depends on the angular frequency Q. For simplicity, we only 
focus on the case of E — injection. If the total effective 
magnetic flux is absent, i.e., On = 0, the transmission and 
reflection are symmetric due to the fact Tr — Ti — T(oj) with 



4*W 



T{u) 



[g* + tHu>-2J) 2 \[g* + tHu + J) 2 ] 



(18) 



Then we have A = Tr - Tl = 0, which characterizes the 
asymmetrical feature. For small effective magnetic flux On « 
0, we have 



A = 



4Jtg 2 T(u>) 
* + t 4 (oj + J) 1 



r*n + 0(On)\ 



(19) 



i.e., the current difference A is a linear function of On in the 
vicinity of On = 0. For a typical case g - yJt, we plot A as 
a function of On in Fig. [3 a). In the resonant case to — 0, it 
shows that this gyroscope have the advantage of good linear 
response within a wide range of On- In practice, in the case of 
extremely large rotation frequency £2, a compensate magnetic 
flux can be added to ensure it works in the linear region. For 
some off-resonant cases, e.g., to - -0.87, the ratio of A to 
On at On = is larger than the one of the resonant case. 
The larger ratio implies a higher sensitivity, so it is good for 
measurement. 

As a result, the instantaneous value Q (t) can be obtained 
by measuring A (r) directly. The cumulative rotation angle is 
just an integration of the angular velocity over time, 



»(t) = 2k a(t)dt. 



(20) 



Therefore, such a system can measure the angular velocity and 
orientation precisely as a gyroscope. 

Actually, the most famous modern gyroscope is the laser 
gyroscope based on Sagnac effect OS El (Two waves prop- 
agating along the opposite directions though a closed rotating 
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FIG. 3: (color online) (a) The difference of the transmission coef- 
ficients on two directions A depends linearly on the total effective 
magnetic flux On around On = approximately. Here, g 2 /t = J, 
E = 0, the frequency detuning to is chosen as to = -0.87 (red dashed 
line), to = (black solid line), and to = 0.87 (blue dot line), (b) 
Intensity distribution of the three-site ring-shaped optical lattice gen- 
erated by two Laguerre-Gauss laser beams. 



ring will induce a relative phase difference, and the position 
of the interference fringes depends on the angular velocity). 
Our proposal is an extention of the Sagnac gyroscope since 
the interference effect is also used. However, the difference 
of the two output currents is measured here instead of observ- 
ing the interference pattern directly, which may enhance the 
sensitivity of the results. 

Our asymmetric coherent scattering system may be realized 
by trapping atoms in a ring-shaped optical lattice. The optical 
lattices with periodic boundary condition have been produced 
by the superposition of two Laguerre-Gauss (LG) laser modes 
theoretically and experimentally 125l |260 . The LG beam with 
frequency to, wave vector k, and amplitude E$ propagating 
along the z axis has the form 



LGf(w) = E /p/(r)^V ( ^- fe) , 



fpiir) = (-IY 



2p\ 
^7T(p + \l\\y 



%(x z )e-^ 2 , (21) 



where x = V2r/ro and ro is the beam waist. L' p (-) is the gener- 
alized Laguerre polynomial. To generate a ID optical lattice 
ring with traps, parameters are chosen as p\ — p2 = 0, 
h = 0, h = N, and X\ = 13x2- The intensity distribution of 
the N = 3 case is shown in Fig. Ob). 



V. GREEN'S FUNCTION APPROACH AND GENERAL 
DISCUSSION FOR ASYMMETRIC TRANSMISSION 

In this section, we generally consider in which case the 
asymmetric transmission can happen for a general multi- 
terminal tight-binding system. For an arbitrary single-particle 
tight-binding model with attaching leads, we will derive the 

T qp between 



condition of the asymmetric transmission T pq + 
lead p and q by using the Green's function formalism intro- 
duced in Ref. B24I1 and references therein. 

We first briefly review Green's function method for multi- 
terminal system. The retarded Green's function of the system 
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as a function of input energy E - —It cos k is 

1 



G R 



E-Hr-T, 



(22) 



C — Pleads 



where He is the Hamiltonian of the central system, and "Lieads 
as the total self-energy denotes the contribution of the half- 
infinite leads, l,i e ads = Upeiiead) £ p . On the single particle basis 

{a\ |0>, a\ |0>, ■ ■ • , a* 10)}, both ff c and S p are N x N matrix. 
Especially, each S p has only one non-zero matrix element at 
the cross point of row p and column p, 



v ' >p,p t 



(23) 



which is because the pth lead is only attached to the pth site of 

' + H.c). Then the advanced 
i [l p - Sj] are 

obtained sequentially. The transmission spectrum from lead p 
to lead q(q + p) is given by 



the central system via (-ga p b l 
Green's function G A 



G R ' and T p matrix r ; , 



T pq =Tr[T q G R T p G A ]. 



(24) 



Applying it to the ring system mentioned in Sec. II-IV, one 
can see that the Green's function formalism gives the same 
exact solution as the one form Bethe Ansatz method. The 
equivalence of the two approaches is proved in the appendix. 

In general, the Green's functions G R (Eq. (l22l ) and G R| are 
complex. From Eq. d24l >. the transmission coefficient T pq is 



- pq 



= Tr [r, (ReG s + /ImG R ) T p (ReG s1 + /ImG Rt )] 
= (T q ) qq Re(G R ) qp (T p ) pp R e (G m ) pq 

-(T q ) qq lm(G R ) qp (r p ) pp lm(G R \ q 
= {T p ) pp (T q ) qq [Re(G% + Im(G%] 

= ^ p)pp(y q)qq \(G ) qp \ , 



where we have use the identity 



Re(G R ) 
Im(G*) 



qp 



Re(G st ) w 
-lm(G m ) 



(25) 



(26) 



together with the property that the F matrix has only one non- 
zero matrix element. Similarly, for the inverse transport, 



Tap — (r P ) PP (r q ) qq \(G R ) 



)pq | 



Therefore, once 



|(G%| = \(G% 



(27) 



(28) 



we have a symmetric scattering T pq = T qp . When this con- 
dition is not satisfied, i.e., |(G S ) 9/ ,| + |(G S ) P9 |, an asymmetric 
scattering is obtained. Eq. ( 122b tells us that the condition is 
not only related to the intrinsic symmetry of the central Hamil- 
tonian H, but also depends on the features of the leads. For 
example, if we only attach two leads on the same 3-site cen- 
tral ring investigated in Sec. Ill, the asymmetric transmission 



phenomena will disappear. This result is in agreement with 
the viewpoint in Ref. [20] that the simplest example must in- 
volve at least three channels. 

Since Hq is Hermitian and 'Luads is diagonal, the term 



Heff - E - H c - Pleads 

can be written as a matrix as 



(29) 



H *ff ~ 



if 

"12 

d* 
"13 



IN 



d\ 2 

d* 

"23 



d 2N 



d\3 
d 2 3 



CN-1 



dlN 



dN-UN 
C N 



(30) 



The diagonal matrix elements are noted by c, while the off- 
diagonal ones are denoted by d. The inverse matrix of H e ff 
has the form 



G R 



D 



Det(H eff ) 



(31) 



where D = Ad}(H e jf) is the adjugate matrix of H e ff, 
which is the transpose of the matrix of cofactors. Since the 
determinant Det(// e yy) is a constant, our problem is reduced 



to the condition of the identity \D, 



•n>\ 



D 



Since a matrix 

and its transpose have the same determinant, we notice that 
D qp and D pq have the form 

Dqp = F\C\, ■ ■ •, Cn\ d\2, ■ ■ ■, dM-\,N\ d* l2 , • ■ ^iv-^Ar) i 

D pq = F\cw; Cn\ d* n , • • •, d* NlN ;d\2, • • ■, ^-i.a') ■ (32) 

Here, F is a polynomial function constructed by c and d. Gen- 
erally speaking, if all of the c and d are non-zero random com- 
plex numbers, \D qp \ is probably not equal to \D pq \. However, 
under some specific conditions, the asymmetric transmission 
will disappear. Now we focus on the resonant systems (there 
is no chemical potential difference between the leads and the 
central lattice) and take some typical cases as an illustration. 

(1) Obviously, when there is no external magnetic flux 
breaks the time reversal symmetry, we have all the d = d*, 
then D qp = D pq , the condition Eq. ( l28T i is automatically satis- 
fied. 



(2) When all the c are real, i.e., c 



D 



D* 



the 

condition Eq. ( f28T > is also satisfied. Because of Eq. ( 1231 ) and 
(TTOl) . this case requires that the incident momentum k is either 
or 7r, which corresponds to injecting an electron at the top or 
bottom of the Bloch energy band. Thus it is not very easy to 
be realized in practice. 

(3) A realizable but not trivial case: The symmetric scatter- 
ing can be obtained even though the magnetic field is present 
if the geometry of the lattice is bipartite and electrons are in- 
jected at E = 0. A bipartite lattice is a lattice that can be 
divided into A and B sublattices, such that a site in A is con- 
nected only with sites in B and vice versa. At this time, the 
function F is a polynomial contains only either even order of c 
or odd order of c. In other words, if one term of F is a product 
of even number of c and some d, no term of F contains the 
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—D; q ifF 

D* pq if F is odd 



product of odd number of c, thus F is called even with respect 
to c. The situation is similar when F is odd with respect to 
c. The momentum for E — is k — n/2, then all the c are 
pure imaginary, c = -c*. Therefore, D qp = (-1)' 

is even with respect to c, while D qp = (-1)' 
with respect to c. For both cases, we have |(G R ) W | = \(G R ) pq ^ 
which represents a symmetric scattering. 

(4) For the bipartite central lattice, T qp equals to T pq even 
though the input energy E + if one of the following condi- 
tions is satisfied, (i) N is odd, there is no lead connected to the 
sublattice with (N - 1) /2 sites, or there are no more than two 
leads connected to the sublattice with (N + 1) /2 sites, (ii) N 
is even, there are no more than one lead connected to either 
sublattice A or sublattice B. 



VI. CONCLUSION 

In this paper, we take the three-site tight-binding ring as an 
explicit example to study the scattering problem of some co- 
herent nanostructure with or without time reversal symmetry 
by virtue of exact solutions. It is discovered that, induced by 
the threaded or effective magnetic flux, the asymmetric trans- 
mission can happen in such a system. With this interesting 
property, on the one hand, such a system behaves as a quantum 
diode for a charged particle due to the pure quantum interfer- 
ence effect. On the other hand, it can also be served as a gy- 
roscope for either a charged or neutral particle when a global 
rotation providing an effective magnetic flux. We have shown 
that the difference of the output currents is linear proportional 
to the rotational angular velocity approximately. The observ- 
able effects presented in this paper are hopeful to be realized 
experimentally in the quantum dot system and the ultra-cold 
atom optical lattice system. We also try to generalize the re- 
sults from the three-site system to the general ones. Accord- 
ing to the Green's function approach B24I1 . it is demonstrated 
that the intrinsic symmetries of the whole system, including 
the geometry of lattice and leads, determine the symmetry of 
scattering coefficients together. Protected by the bipartite con- 
figuration, a symmetry scattering may be obtained even if the 
time reversal symmetry is broken for the central system. 



The Hamiltonian of a general central system with N sites 
reads 



H C = ^ (//(•);.,• (l)a r 

jj'=l 



(A.l) 



We suppose that there is a set of j defined as {lead} = 
0'i»j2. Jm), where 2 < M < N. If e {lead}, a half 
infinite tight-binding lead is attached to the jth site. The total 
Hamiltonian of leads is 



je{lead| 



H L = - Sa)bf +tY J b^ j] bf +i +H.c. 



(A.2) 



According to the Bethe Ansatz method, the scattering state is 
supposed to be 

l<A> = X> ' )fl } |0>+ Z Tj^Ulb^lO), (A3) 



where 



&,Sf) -B(j)e ikl ' -6 „,<■ (A A) 



The pth lead has been chosen as the input lead. From the 
Schrodinger equation 



(H C +H L )\M = E\M, 



(A.5) 



i// a (j) and ijrr (j) should satisfy 



2 (H c )jj i/f a (/) - gtfjttad}^ 1 (1) = Eiff a (j) (A.6) 



for any j, and 



-t (f + 1) + ^ (f - 1)] = E^ ] (/) (A.7) 
for j € {lead}. Eq. ( IA.4I ) and ( lA.7b give the dispersion relation 
£ = -2fcos/t, (A.8) 
together with the equation for j e {lead}, 
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APPENDIX: EQUIVALENCE BETWEEN BETHE ANSATZ 
METHOD AND GREEN'S FUNCTION METHOD 



In this appendix, we prove that the Bethe Ansatz method 
and Green's function method are equivalent for determining 
the transmission coefficients. 



^ a (j) = B(j) + dj, p . (A.9) 
For future convenience, we define a column vector X, 
X - (Xi X2 ■ ■■ Xn) , 

Xj = (l-^fleadjf^aOlj + ^iead^C/). (A- 10) 

Eq. iA.6i can be simplified as 



Z 



E6jj - (H c )jj, -\-—e' 



(Hc)j, p -Sj, P \ E + -< 



(A.ll) 
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which is a system of linear equations in the variables Xj, In 
the matrix form, Eq. jA. lib becomes 

H eff X=W, (A. 12) 

where 

Heff - E - H -Pleads- (A. 13) 

The matrix elements of H e /f, W, and £/ ea dj are 

(Heff) jf = E8JJ -{H c )jj - \- 8 jj%j6 jA ^, 



and 



(so-s;) 



(A. 18) 



is a real number. So that 



Det (H® f ) = Del = ir A pq . (A. 19) 



A M is the (p, q)th algebraic cofactor of ff e //, which is actually 
{—\) p+q times the determinant of the submatrix obtained by 
removing pth row and qlh column from H e ff. Therefore, the 
transmission coefficient from lead p to lead q is 



(Eleads) j,f - -— e' 5;,/5/{lead) - ^0<J>,/^',(lead), (A. 14) 

where Eq = - (g 1 It) e' k . The solution of Xj can be given by 
the Cramer's rule Il27ll . Without loss of generality, we focus 
on the solution of X q with q e {lead} but q + p, 



X q = B(q) = 



Det(tf e// )' 



(A. 15) 



Here, the matrix elements of /r are 

eff 

( H eff)jj - ("-"),., + 6 M W r (A. 16) 

The determinant of ff^// equals to the one of another matrix 
/T,l where 

- (i - */, ? ) + «^r , (A. l?) 



[ pq 



\B(q)\ 2 



TnA 



O^pq 



Det(ff e// ) 



(A.20) 



The definition of the inverse of a matrix tells us that Eq. ( IA.6I 1 
is actually 



^= r o|(^//)J = r o|( G \| • (A.21) 
On the other hand, from Green's function method Eq. d25l ). 



Tpq - 



(A.22) 



where (J p) PP = (E 9 ) f/9 = To. Then Eq. ( 1A.22I ) and ( 1A.21I ) are 
same, the Bethe Ansatz method is equivalent to the Green's 
function method for determining the transmission coefficients. 
Moreover, the Bethe Ansatz method can also gives the exact 
wave function. 
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